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' Abstract We study the decay dynamics of an inter- 
acting Bose-Einstein condensate in the presence of a 
metastable trapping potential from which the conden- 
. sate can escape via tunneling through finite barriers. The 
' time-dependent decay process is reproduced by means 
, of the instantaneous decay rates of the condensate at a 
' given population of the quasi-bound state, which are cal- 
culated with the method of complex scaling. Both for the 
case of a double barrier potential as well as for the case 
of a tilted periodic potential, we find pronounced devi- 
ations from a monoexponential decay behaviour which 
would generally be expected in the absence of the atom- 
' atom interaction. 



1 Introduction 

With the advent of optical lattices [Tl[2l[3] and "atom 
chips" [HIS], it became possible to probe the transport 
properties of a Bose-Einstein condensate in the meso- 
scopic regime. The unprecedented degree of experimen- 
tal control in these systems led to the observation of 
Bloch oscillations [6l|7] , the guided and free propagation 
of condensates through waveguide structures [HUS], the 
transport of condensates with "optical tweezers" [TU] , as 
well as the realization of Josephson junctions [TT] and 
matter- wave interferometry |12| . to mention just a few 
examples. Those experiments typically involve rather small 
trapping potentials, with length scales that can be of 
the order of a few microns. In such geometries, decay 
mechanisms of the condensate become a relevant issue. 
On the one hand, the condensed state is, at finite atom 
densities, subject to depletion, which is caused by the 
interaction with the thermal cloud and by three-body 
collisions. On the other hand, the condensate can escape 
from the trapping potential by tunneling though its bar- 
riers, if the chemical potential of the condensed atoms 
exceeds the background potential in the free space out- 
side the trap. In that case, the self-consistent mean-field 



state of the condensate is no longer bound, but rather 
corresponds to a metastable "resonance" state, in a sim- 
ilar way as, e.g., doubly excited electronic states in the 
helium atom [13j . 

From the theoretical point of view, various methods 
were used [T4l[T5l[T6lfT7ra fT9] to tackle the problem of 
how to treat "resonances" , i.e. stationary states that de- 
scribe the escape of population from an open confine- 
ment potential, in the context of Bose-Einstein conden- 
sates. For linear systems, it is well known that this task is 
most conveniently accomplished by applying the method 
of "complex scaling" (or "complex rotation") [20ll2ip 22 [ 
I23j . This technique essentially amounts to the complex 
dilations r i-^ re'^ and — iV i— > — iVe^*^ of the posi- 
tion and momentum operators in the Hamiltonian that 
describes the quantum system under study. This trans- 
formation leads to a nonhermitean Hamiltonian with 
a complex eigenvalue spectrum the continuous part of 
which is rotated to the lower half of the complex energy 
plane. Resonances, i.e. decaying states with eigenvalues 
corresponding to poles of the resolvent below the real en- 
ergy axis, are thereby uncovered and can be calculated 
using standard diagonalization techniques for complex 
matrices. This approach is essentially exact, in the sense 
that no a priori approximations are introduced in the 
complex dilation procedure. Positions and widths of res- 
onances can therefore be calculated with high precision 
by means of the complex scaling procedure 22^23j . 

The generalization of this approach to Bose-Einstein 
condensates was recently accomplished in our previous 
studies jT8lfT9| , where we apply the complex scaling trans- 
formation to the nonlinear Gross-Pitaevskii equation that 
describes the mean-field dynamics of the condensate. In 
contrast to an alternative approach proposed by Moi- 
seyev and Cederbaum |17j, we explicitly take into ac- 
count the complex nature of the wavefunction of the res- 
onance state, which leads to a considerable complication 
of the problem due to the resulting nonanalyticity of the 
interaction term in the Gross-Pitacvskii equation. We 
showed in Ref. \18\ how this complication can be tackled 
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and how quasi-bound resonance states of the condensate 
can be calculated by means of a real-time propagation 
approach based on the complex-scaled Gross-Pitaevskii 
Hamiltonian. 

In this paper, we apply this method in order to cal- 
culate specific time-dependent decay processes of the 
condensate. Instead of a direct numerical integration of 
the time-dependent Gross-Pitaevskii equation (which is 
rather consumptive in CPU time for the evolution time 
scales under consideration), we compute, with complex 
scaling, the decay rates of the quasi-bound state of the 
condensate at various values for the effective interaction 
strength (which would be proportional to the number of 
atoms that are populating this quasi-bound state at a 
given instance of time). Then we integrate, on the basis 
of this information, a simple rate equation that directly 
describes the decay of quasi-bound population. This ap- 
proach is comparatively efficient, avoids the introduction 
of artificial complex potentials at the grid boundaries 
in order to absorb the outgoing population (see, e.g., 
Ref. [Ml), and provides physical insight that could be 
used to control the decay process in a similar way as 
for macroscopic tunneling of condensates in double well 
potentials (e.g., ^24,25J). 

The paper is organized as follows: In Section[21 we es- 
tablish the general relation between resonance states of 
the stationary Gross-Pitaevskii equation and the actual 
time-dependent decay process of the condensate. We fur- 
thermore discuss how such resonance states can be calcu- 
lated by the method of complex scaling, as was described 
in more detail in Ref. [18] . Section[3]contains the numeri- 
cal results that we obtain for two paradigmatic examples 
of metastable confinement configurations: a harmonic 
trapping potential with Gaussian envelopes, and a tilted 
periodic lattice. We calculate the time-dependent decay 
of the condensate, which is, in both cases, characterized 
by a pronounced nonexponential nature. 

2 The nonlinear complex scaling approach 

We consider a Bose-Einstein condensate that is confined 
within a cylindrical matter-wave guide with transverse 
frequency uj±_ and evolves in the presence of a longitudi- 
nal potential V{x). In the "ID mean-field regime" [26| 
(where the confinement is strong enough to inhibit trans- 
verse excitations within the waveguide, but not as strong 
as to enter the Tonks- Girardeau regime [21]), the dynam- 
ics of the condensate is described by the one-dimensional 
time-dependent Gross-Pitaevskii equation 

d / 1 \ 

i-V'(x,t) = f --_ + y(a:) +5o|^(a;,i)n ^i.x,t) 

(1) 

where x denotes the coordinate along the waveguide. The 
longitudinal potential V{x) is assumed to provide a local 
harmonic confinement with trapping frequency lo\\^ from 
which the condensate can escape via tunneling through 



finite barriers. Dimensionless variables, defined by set- 
ting fi = m = =1, are used throughout this paper. 
This means that we express all length scales (including 
the coordinate x) in units of a|| = ^?i/(mw||), all energy 
scales in units of and all time scales in units of 
The effective one-dimensional interaction strength is, in 
these units, given by go — SosLJj^ where denotes the 
s-wave scattering length of the atoms . 

For the description of our theoretical approach, we 
specifically focus in the following on the double barrier 
potential 

\/(a;) = ia;2exp(-ax2) (2) 

with a = 0.1, which could be experimentally realized, 
e.g., with red- and blue-detuned laser beams that are 
tightly focused onto the waveguide. Obviously, V{x) does 
not exhibit any bound state, and the eigenspectrum of 
the linear (noninteracting) Hamilton operator is fully 
continuous. There exist, however, quasi-bound states which 
are localized in the well around a; = and which give rise 
to resonances in the energy spectrum (corresponding to 
complex poles of the scattering matrix). In the case of 
noninteracting atoms {go = 0), such resonance states are 
described by time-dependent wavefunctions of the form 
ipix^t) — ip{x) exp{~iEt), where ip{x) satisfies the sta- 
tionary Schrodinger equation for the complex eigenvalue 
E = fi — ir/2 and exhibits outgoing (Siegert) boundary 
conditions [28] tpi^) — ^ V'o exp(i/c|x|), with Re(fc) > 
for X ±oo. This latter property expresses the fact that 
the wavefunction of the decaying state is characterized 
by a finite current of atoms that propagate away from 
the well. 

As a consequence, the atomic population inside the 
well decays exponentially according to oc exp(— I^t) if 
the system is initially prepared in the energetically low- 
est resonance state. It is quite obvious that this is no 
longer true in the nonlinear case of interacting atoms 
(ffo 7^ 0). There, the tunnel coupling through the bar- 
riers explicitly depends, via the nonlinear term in the 
Gross-Pitaevskii equation, on the local atomic density, 
which in turn induces a temporal variation of the de- 
cay rate F. We therefore naturally obtain, as was also 
pointed out in Ref. [29] . a nonexponential decay of the 
atomic density, the reproduction of which is the central 
aim of this paper. 

Despite this complication, a description of the decay 
process of the interacting condensate in terms of instan- 
taneous quasi-bound states can nevertheless be justified 
if the rate F characterizing the temporal variation of the 
density inside the well is rather small compared to the 
chemical potential (which should generally be the case 
if the condensate escapes via tunneling through finite 
barriers). In such a quasi-stationary situation, we can 
employ an adiabatic ansatz where the condensate is as- 
sumed to remain always in the energetically lowest (and 
most stable) resonance state associated with a given in- 
stantaneous density \'ip{x,t)\'^ . This resonance state is 
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formally defined, together with its associated complex 
eigenvalue Eg = fj,g — irg/2, by the self-consistent solu- 
tion of the nonlinear stationary equation 

Hi,Pg)i;gix) = Egt^gix) , (3) 

with 

HW^-l^ + V{x)+gmx)\' . (4) 
ipg{x) is normalized according to the condition 

/oo 
\i;g{x)\^w{x)dx = l. (5) 
-oo 

Here u'(x) represents a weight function that measures 
the population inside the well. For the double barrier 
potential a natural choice for this weight function 
would be w{x) = 9{a — x)6{x + a) where a = 
corresponds to the maximum of the potential and 0{x) 
denotes the Heavyside step function. 

In addition, tpg{x) should also satisfy outgoing bound- 
ary conditions, which would imply an asymptotic be- 
haviour of the form ipgi^x) (x exp[z J k{x')dx'] for large 
a; ^ oo (and a similar one for x —^ —oo) where the 
spatial dependence of the effective wavenumber k{x) ac- 
counts for the smooth variation of the self-consistent po- 
tential in Eq. ([4]) . It was pointed out in Ref . TS' that this 
condition can only be fulfilled in an approximate way up 
to a given maximum spatial distance Xc at which the in- 
teraction energy g\i(jg{x)\'^ starts to exceed the chemical 
potential /i. For x > Xc, the self-consistent quasi-bound 
state would formally encounter a singularity, which re- 
flects the fact that explicit time-dependence is expected 
beyond that critical distance. 

On the basis of these resonance states, we can now 
formulate the adiabatic ansatz for ip(x,t) as 

tp{x,t) ^ ^/N^^/jg(^t){x)exp (^-i J Eg(t')dt'^ ■ (6) 

Here, the effective time-dependent interaction strength 
is given by g{t) = goN{t) where N{t) denotes the time- 
dependent population inside the well, defined by 

/oo 
\'4!{x,t)\^w{x)dx. (7) 
-oo 

Using the normalization condition ([5]) of the resonance 
state and taking into account the fact that its eigenvalue 
Eg = fig — irg/2 IS complcx, one can straightforwardly 
derive the implicit expression 

N{t)^Noexp(^^ rg^t,^dt'Y (8) 

for the quasi-bound population of the condensate, which 
can also be formulated in terms of the ordinary differen- 
tial equation 

^ = -rg^t)Nit) , (9) 



with the initial condition A^(0) — Nq. The time-dependent 
decay process of the condensate can therefore be entirely 
reproduced with comparatively little numerical effort if 
the decay rates -^'^(t) of the instantaneous quasi-bound 
states 4'g(t) (x) are known. 

As was described in detail in Ref. 'JSj , the calculation 
of the decay rates can be achieved by the method of 
complex scaling. This technique essentially amounts to 
the application of the nonunitary mapping 

^{x) ^ ■4}^^\x) = Rs^j{x) = e'^l'^i){xe'^) (10) 

to the wavefunction ip{x), which corresponds to the com- 
plex dilation x i— > xe^^ of the position operator. Applying 
this transformation to the linear stationary Schrodinger 
equation Hqi/j — Eip — with Hq being defined through 
Eq. ^ via Hq = H{ij} = 0) — yields the complex sta- 
tionary equation 

H^^\(<>\x) = Eip'^'^\x) (11) 
with the complex scaled Hamiltonian 

ff^ ^ ReHoRe' = -\^~^'' ^ + ^(^e^') • (12) 

The spectral properties of this nonhcrmitcan Hamil- 
tonian are widely discussed in the literature on com- 
plex scaling [20 l [2T 1 l22 l [23] : While bound states of the 
original Hamiltonian Hq (which are absent in our par- 
ticular case) remain bound after the complex dilation 
(as long as \9\ < 7r/4), the continuum states are "ro- 
tated" in the complex energy plane, in such a way that 
their eigenvalues are located along the axis E = ee~^'^ 
with real positive e. This rotation uncovers the spec- 
tral resonances of the system, which correspond to the 
poles of the analytical continuation of the Green func- 
tion G = {E — Hq + iS)~^ to the lower half part of 
the complex energy plane. Those resonances turn into 
discrete complex eigenvalues En — fin — iEn/^ under 
complex dilation, and are represented by normalizable 
eigenfunctions that can be straightforwardly calculated 
by diagonalizing in any numerical basis. 

The generalization of this approach to the nonlin- 
ear case would be comparatively straightforward if the 
replacement |?/'(x)p ['ijj(x)Y in the nonlinear Gross- 
Pitaevskii Hamiltonian ^ could be justified. For this 
particular case, the implementation of the complex scal- 
ing technique was explained in detail in Ref. [T7]. In 
reality, however, the wavefunction of the resonance state 
is intrinsically complex due to the outgoing boundary 
conditions (i.e, ^j{x) cx exp(ifc|a;|) for — > oo), and the 
resulting nonanalyticity in the Hamiltonian ^ intro- 
duces a major complication of the problem. Formally, a 
second analytic wavefunction ip needs to be introduced, 
which coincides with the complex conjugate of on the 
real axis, i.e. 

il^ix) = '4>*{x) for real x, (13) 
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and which is independently transformed under the nonuni- 
tary dilation operator Rg^ i.e. 

i^ix) ^ i^^'^\x) EE Re^(x) = e*^/2^(xe^^) . (14) 

The analytic continuation of the stationary Gross-Pitaevskii 
equation to the complex domain yields then 

H^^\il))il;^'^\x) = Eil;^'^\x) (15) 

where the complex scaled nonlinear Hamiltonian is given 

by 

H^'^Xi^) = H^^^ + ge-''^^^'^\x)xl:^'^\x) . (16) 

The lowest resonance state of the condensate can 
be calculated by a real-time propagation approach [T5] . 
i.e. by numerically propagating -0^*^ under the time- 
dependent Gross-Pitaevskii equation 

i^^^^\x) ^ H^'\4,,)^^^\x) (17) 

in the complex scaled system, where r represents a ficti- 
tious numerical "time" parameter (which is unrelated to 
the physical time evolution in the actual decay process). 
In practice, V'r {x) is expanded on a spatial grid, and 
an implicit finite-difference scheme is employed to carry 
out the mapping ■f/^i^-' ^ i'^+sr small time steps St. If 

(9) 

ijjr is renormalized after each propagation step in order 
to satisfy the condition (O, the integration of Eq. (flT)) 
necessarily converges to the most stable resonance state 
of the complex scaled Hamiltonian, which corresponds 
to the quasi-bound state with the smallest decay rate 

m- 

The major numerical difficulty in this approach lies 
in the evaluation of the nonlinear term in the complex 

scaled Hamiltonian p6p . Indeed, V (^) is not identi- 
cal to [tp^^^x)]* , the complex conjugate of ^'''^■'(2;), and 
needs to be evaluated according to the relation 

lp'-'\x) ^ Re (rI^^ (x) (18) 

which requires explicit back- and forward-rotations of 
the complex scaled condensate wavefunction. In prac- 
tice, these rotations (which are also used to evaluate the 
normalization condition ([5])) are numerically performed 
by mapping the grid representation of the wavefunction 
into a nonorthogonal set of analytic Gaussian orbitals 
01, (x) that are centred around different positions along 
the grid, and by using a transformation matrix that 
contains the overlap integrals / 4)^{xe'^^)(t>^i {x)dx as ele- 
ments (see Ref. [18] for more details) . Such an operation, 
however, is known to be potentially unstable j31] and re- 
quires great care in the numerical implementation. It is 
therefore not obvious to which extent unlimited preci- 
sion in the decay rates of the quasi-bound states can be 
achieved within this nonlinear complex scaling approach. 



3 Calculation of time-dependent decay processes 

3. 1 Double harrier potential 

Despite this latter complication, we find that the chemi- 
cal potentials and decay rates of the self-consistent quasi- 
bound states of the double barrier potential ^ can be 
calculated in this way with rather good accuracy, even 
in case of very strong nonlinearities where the resonance 
level lies close the barrier height of V{x). This was ex- 
plicitly verified in Ref. [18j by comparing the resulting 
values for /ig and Fg with the ones that are obtained from 
an alternative approach, which was based on the real- 
time propagation of the original (i.e., unsealed) Gross- 
Pitaevskii equation in the presence of absorbing bound- 
aries. Good agreement was generally found between the 
two approaches [TB]. 

The chemical potentials and decay rates of the low- 
est resonance state are plotted in Fig. [T] as a function 
of the effective interaction strength g. Quite intuitively 
/ig increases with increasing g due to the presence of 
the mean-field interaction energy in the nonlinear Gross- 
Pitaevskii Hamiltonian. This increase of the chemical 
potential results in a dramatic enhancement of the de- 
cay rate Fg, which can be explained by the fact that 
the effective imaginary action integral that semiclassi- 
cally determines the tunneling rate through the barri- 
ers is appreciably reduced with increasing energy. It was 
pointed out in Refs. [M l flB l fTS] that an attractive inter- 
action between the atoms leads to a stabilization of the 
resonance state, i.e., to a reduction of the chemical po- 
tential to values below /i = where the associated decay 
rate would vanish. For the double barrier potential un- 
der consideration, this stabilization process would occur 
at g ~ —1.1. 

With this information, we can now quantitatively re- 
produce the time-dependent decay process of the con- 
densate by means of the integration of the rate equation 
([9]). For this purpose, we use the values of the decay rates 
that are calculated with the complex scaling method at 
the equidistant interaction strengths g — 0,0.1,0.2..., 
and employ a cubic interpolation to obtain intermediate 
values of F. As initial value of the effective interaction 
strength, we consider g(t = 0) = 4, where the chemi- 
cal potential of the quasi-bound state lies already rather 
close to the barrier height of the potential. In the spe- 
cific case of a condensate of ^^Rb atoms that encounters 
the longitudinal and transverse confinement frequencies 

— u}±^ ~ 2tt X 10^^ Hz, this would imply that about 
iV(0) = 4/(70 — 100 atoms are initially localized in the 
single well [32] . 

The result of the integration is displayed in the two 
lower panels of Fig. [1] We see a clearly nonexponential 
decay of the bound population, which refiects the fact 
that the decay rate decreases with descreasing interac- 
tion strength g. As a consequence, a rather large num- 
ber of atoms leave the trap during the first 1000 units 
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Fig. 1 Decay of a Bose-Einstein condensate in the double 
barrier potential. The upper panel shows the potential ^ to- 
gether with the chemical potential of the initial quasi-bound 
state of the condensate (thick horizontal line), which corre- 
sponds to the initial value g{t = 0) = goN{t = 0) = 4 of the 
effective interaction strength. During the decay process, the 
chemical potential decreases with time, due to the reduced 
interaction energy, and approaches the level of the noninter- 
acting quasi-bound state (thin horizontal line in the upper 
panel). This decrease of /i is accompanied by a strong reduc- 
tion of the decay rate F, which can be seen in the two middle 
panels where /i and F are plotted as a function of g. As a 
consequence, a pronouncedly nonexponential decay for the 
quasi-bound population P{t) = N{t)/N{t = 0) is obtained 
from integrating the equation dN/dt = —r{N)N, which is 
displayed in the two lower panels (where P is plotted vs. t 
on a linear and logarithmic scale, respectively). 



of the evolution time, while the remaining part of the 
condensate becomes stabilized and decays, for asymp- 
totically large times, with the rate Iq ~ 2 x 10~^ of 
the noninteracting quasi-bound state. In the example of 
a ^^Rb condensate in a confinement with trapping fre- 
quency L^ii = 27r X 10^ Hz, the above characteristic time 
scale of the nonexponential behaviour would correspond 
to t ~ 100 ms which is of experimental relevance. 

3. 2 Tilted periodic potential 

Nonexponential features can also be observed in the pres- 
ence of relatively small interaction strengths, namely if 
the confinement potential permits the possibility of res- 
onant tunneling. This is, for instance, the case for the 



tilted periodic potential 

V{x) ^sm^{x/2) + Fx (19) 

that is experimentally realized with optical lattices [lll2l 
13] or, within the atom chip context, by means of periodic 
sequences of microfabricated wires [3 3) . In this nonlinear 
Wannier-Stark system, the possibility of resonant tun- 
neling arises if the local ground state of one of the wells 
is nearly degenerate with the first excited state of the 
adjacent well. As was pointed out in Ref. [55], this near- 
degeneracy would give rise to a significant enhancement 
of the condensate's decay rate. 

As in our previous study [12], we assume that the 
condensate is initially confined within one single well of 
the lattice. The method of complex scaling can again be 
used to calculate the chemical potential and decay rate 
of the self-consistent quasi-bound state in the presence 
of the interaction, even though the tilted potential ([TO)) 
leads to an asymptotic spatial behaviour of the contin- 
uum states that is substantially different from the previ- 
ous double-barrier problem. As was described in detail in 
Ref. [in] , additional complications arise in this potential 
(such as the existence of many different self-consistent 
resonance states with identical decay rates) and tech- 
nical modifications need to be implemented in order to 
achieve good convergence of the real-time propagation 
method. In analogy with the double barrier potential, 
the weight function that characterizes the bound pop- 
ulation inside the well according to Eq. ^ is given by 
'w(x) = 9{'K — x)9{x + tt). 

The resonance-enhanced decay process of the con- 
densate in this tilted lattice is displayed in Fig. [21 for the 
tilt strength F — 0.1412 at which the level of the non- 
interacting local ground state in each well lies slightly 
above the level of the first excited state in the adjacent 
well on the left-hand side. As initial value for the nonlin- 
earity, we consider g(t = 0) = —0.25, which would cor- 
respond to an attractive interaction between the atoms 
(which could be realized, e.g., by using condensates with 
''Li atoms or by applying Feshbach tuning techniques 
[34]). This attractive nonlinearity lowers the chemical 
potential of the self-consistent quasi-bound state, in such 
a way that is becomes shifted below the level of the first 
excited (noninteracting) state in the adjacent well. Dur- 
ing the time evolution, the loss of population leads to 
an increase of the chemical potential, which, at about 
t ~ 2500 time units, approaches the resonance. At that 
point, the decay rate becomes drastically increased re- 
sulting in a significant escape of atoms from the well. As 
a consequence, the chemical potential of the quasi-bound 
state quickly moves out of resonance, and the decay rate 
becomes again reduced. 

The time evolution of the condensate was again cal- 
culated by integrating the rate equation ([5]), using de- 
cay rates that were calculated with the complex scal- 
ing method at the equidistant values <? = 0, —0.0025, 
—0.005, ... of the interaction strength. The accuracy with 
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Fig. 2 Decay of a Bose-Einstein condensate in tlie tilted 
periodic potential (|19|l with F — 0.1412. The atomic cloud 
is assumed to be entirely localized within a single well of the 
lattice, and decays via tunneling through the barrier on the 
left-hand side of the well. This decay process can be sub- 
stantially accelerated in the case of resonant tunneling, i.e., 
if the chemical potential of the condensate matches the en- 
ergy of the first excited (noninteracting) state of the adja- 
cent well (dashed horizontal lines; the thin solid lines mark 
the noninteracting ground levels of the wells). In contrast 
to the case of the double barrier potential, we consider here 
a relatively weak and attractive initial interaction strength, 
g{0) = —0.25, which leads to an initial chemical potential 
that lies slightly below the level of this excited state (thick 
horizontal line in the upper panel). The time evolution of 
the quasi-bound population P{t) = N{t)/N{0), shown in the 
lower panels, clearly displays the characteristic signature of 
an intermediate resonant tunneling process. Note that the 
wiggles in the decay rates at large \g\ (middle right panel), 
which arise from numerical inaccuracies, do not leave signif- 
icant traces in the time evolution of P{t). 



which these decay rates could be calculated was not as 
good as in the case of the double barrier potential [18 , 
which is clearly reflected by the appearance of wiggles 
on the left-hand side in the middle right panel of Fig. O 
These artificial fluctuations of F, however, leave no sig- 
nificant traces in the actual time-dependent decay pro- 
cess, as is clearly seen in the lower two panels of Fig. [21 
This intermediate "burst" of atoms should be readily 
observable within existing experimental setups based on 
optical lattices or atom chips. By imposing a rather weak 
transverse confinement, the one-dimensional interaction 
strength (given by go = 2ashuj± [27]) can be appreciable 



reduced, what should allow one to suppress effects be- 
yond the mean-field description of the condensate, and 
additional longitudinal potentials can, as in the exper- 
iment on tunneling well potential [35] . be employed to 
prepare the condensate in one single well of the lattice. 
Obviously, the attractive interaction between the atoms 
is not a necessary condition for this nonexponential de- 
cay phenomenon: Indeed, the same effect could be in- 
duced with a repulsively interacting species at slightly 
weaker tilt strengths F, where the noninteracting ground 
state of the well lies slightly below the first excited level 
of the adjacent well. In both cases of attractive and re- 
pulsive interaction, the intermediate enhancement of the 
decay rate should clearly manifest in the time-of-flight 
image of the condensate after the decay process, which 
would display a pronounced peak due to the effect of 
resonant tunneling. 



4 Conclusion 

In summary, we studied the time-dependent decay of 
Bose-Einstein condensates in mesoscopic trapping po- 
tentials that permit escape by tunneling through finite 
barriers. The decay process of the condensate was repro- 
duced by integrating a simple rate equation for the quasi- 
bound population, using instantaneous decay rates that 
were computed by means of the method of complex scal- 
ing. This approach is rather efficient as compared to a di- 
rect numerical integration of the time-dependent Gross- 
Pitaevskii equation, and provides additional insight into 
the mechanisms that underly the decay of the conden- 
sate. Though only applied for one-dimensional configu- 
rations, the complex scaling approach can be straightfor- 
wardly generalized to three-dimensional decay problems, 
and might furthermore represent a convenient concep- 
tual framework for treating resonances of the nonlinear 
Gross-Pitaevskii equation from the mathematical point 
of view. 

With this approach, we calculated the decay of a 
Bose-Einstein condensate in a double barrier potential 
P!^[TMT5] as well as in a tilted periodic potential [IHl 
[T^ . For this latter case, we find a strong intermediate 
enhancement of the tunneling rate, which arises due to 
a near-degeneracy with a quasi-bound state in another 
well of the periodic potential. This enhancement leads to 
a pronounced deviation from an exponential behaviour 
of the condensate's escape dynamics, which could be 
controlled by suitable time-dependent variations of the 
tilt field, in a similar way as for macroscopic tunneling 
in double well potentials [24ll25] and for nonlinear reso- 
nant transport through atomic quantum dots 36J. Such 
nonexponential effects should be readily observable in 
present-day state-of-the-art experiments on interacting 
matter waves in mesoscopic trapping potentials [51 1111 
a2ii33l[33[37]. 
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